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Abstract 

We study solutions of the Yang-Baxter equation on a tensor product of an arbitrary finite¬ 
dimensional and an arbitrary infinite-dimensional representations of the rank one symmetry algebra. 
We consider the cases of the Lie algebra 5^2, the modular double (trigonometric deformation) and 
the Sklyanin algebra (elliptic deformation). The solutions are matrices with operator entries. The 
matrix elements are differential operators in the case of si 2 , finite-difference operators with trigono¬ 
metric coefficients in the case of the modular double, or finite-difference operators with coefficients 
constructed out of Jacobi theta functions in the case of the Sklyanin algebra. We find a new factorized 
form of the rational, trigonometric, and elliptic solutions, which drastically simplifies them. We show 
that they are products of several simply organized matrices and obtain for them explicit formulae. 
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1 Introduction 

A quantum integrable system corresponds to each solution of the Yang-Baxter equation 

Ruiu - v) Mi 3 (u) R23{v) = K23(?^) - v) . (1) 

Consequently classification of its solutions is of major interest for mathematical physics. The linear 
operator Mij(n) from Eq. ([T]) that acts on a tensor product of two linear spaces and depends on a 
spectral parameter u € C is traditionally referred to as H- matrix. We prefer to call it Ji- operator, since 
we will extensively work with infinite-dimensional linear spaces. The operators in Eq. © act on a 
tensor product of three linear spaces. Each R-operator acts non-trivially on a pair of spaces denoted by 
its lower indices, and it is extended as an identity operator on the remaining space of the triple. 

It is well known that solutions of the Yang-Baxter equation can be rather intricate [221126] . None 
the less appealing to the Quantum Inverse Scattering Method [l5l[28] one can put forward a reasonable 
conjecture that they are composite objects having internal structure and that they are constructed out 
of elementary blocks. A more refined statement is that the R-operator admits factorization, i.e. it 
is a product of several simpler operators. This observation enabled to construct the general solution 
of the Yang-Baxter equation,Eq. ©, acting on a tensor product of two infinite-dimensional principal 
series representations of the group SL(A^, C) [12] . In the case of rank one algebra this result has 
been carried over to trigonometric and elliptic deformations. The general R-operators for the Faddeev’s 
modular double and the elliptic modular double has been constructed in a factorized from in [7] and m, 
respectively. 

In this note we will deal with finite-dimensional representations. We will prove that R-operators for 
rank 1 algebras acting on a tensor product of an arbitrary finite-dimensional and an arbitrary infinite¬ 
dimensional representations admit factorization as well. These solutions of the Yang-Baxter equation, 
Eq. ©, can be thought of as generalizations of the quantum Lax operator, since the fundamental 
representation in the auxiliary space is substituted by a higher-spin representation in 

Let us consider firstly solutions of Eq. ([T]) that are invariant with respect to the Lie algebra s^ 2 ■ In 
the following sections we will consider as well its trigonometric deformation that is the modular double 
(along with ?7g(sQ)) and its elliptic deformation that is the Sklyanin algebra. 

The commutation relations between s ^2 generators are the following 

[S+,S-] = 2S, [S,S^] = ±S^. (2) 

The symmetry restriction implies commutativity of the R-operator and the co-product of the algebra 
generators 

[Ru{u),Sf + Sf]=0 , []Ri 2 (u), S 1 + S 2 ] = 0. 

The linear spaces the R-operator acts upon are representation spaces of sQ. We will be concerned with 
representations of si 2 that are Verma modules. We realize Verma modules in the space of polynomials 
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C[z]. The generators of st '2 are first order differential operators acting on the space of polynomials and 
depending on a parameter £ € C, which we call spin of the representation, 

S = zd-i , S- = d , S+ = -z‘^d + 2ez . (3) 


At generic i the action of the generators ([3]) on the space ^[z\ is irreducible, and the representation 
is infinite-dimensional. At (half)-integer spins 2£ = re € Z>o the representation is reducible, and a 
(re + l)-dimensional irreducible representation with the basis {1, z, z^, ■ ■ ■ , z”} decouples. We prefer to 
gather the basis monomials in the generating function (z — x)^ which depends on an auxiliry parameter 
X. Expanding the generating function with respect to x we recover all basis vectors. 

An elegant formula for st' 2 -invariant solutions of the Yang-Baxter equation, Eq. m, acting on a 
tensor product of two representations of arbitrary spins s and i has been derived in [T91I77] . 


where J is a “square root” of the Casimir operator: J(J + 1) = (5i -|- 82 )'^; the operator P 12 swaps 
the tensor factors: Pi 2 ‘h( 2 ;i, 2 : 2 ) = ^(z 2 ,zi). The formula ([H) is valid for both finite-dimensional and 
infinite-dimensional representations. The operator J is defined rather formally. Thus the formula ([4]) has 
to be accompanied with a decomposition of the tensor product of two representations into irreducibles, 
which are eigenspaces of J. 

In [23] the universal R-matrix for the Yangian double of 8^2 has been found in a form of a product 
of three power series in generators S, S^, Eq. ([2|). This universal R-matrix taken in the evaluation 
representation is an alternative to Eq. (HD- 

Here we choose another opportunity. We will obtain a number of explicit formulae for solutions of 
the Yang-Baxter equation, Eq. working with the functional realization of representations, Eq. Q. 
Indeed, the R-operator acting on the space of polynomials of two complex variables C[ 2 ;i] (8) C[ 2 ; 2 ] takes 
the form 


Ri2(re|s,£) = P 12 


r(2:2i<92 - 2s) 
T{z 2 id 2 - u- s - t) 


T{zi 2 di + u- s - t) 
T{zi 2 di - 2 s) 


( 5 ) 


where Zij = Zi — Zj. We imply that representations of the form (|3|) specified by spins s and i are realized 
in the spaces Cfzi] and C[ 2 ; 2 ], respectively. The ratio of two gamma functions of operator argument can 
be rewritten as an integral operator by means of the Euler integral of the first kind 


T{zi 2 di -F a) 

r(zi2di + b) 


^(zi,Z2) 


1 

r(6-a) 



daa'^ ^(1 —“ ^^{azi + {l 


a)z2.,Z2) ■ 


Let us note that the R-operator in Eq. (|5|) is factorized. The origin and the meaning of this and 
other similar factorizations has been clarified in m- The equality of R-operators (jH) and ([S|) (up to an 
inessential normalization factor), provided the functional realization of 5 ^ 2 , Eq. Q, is adopted, can be 
checked by a straightforward calculation |12j . 

The solution ([5|) of the Yang-Baxter equation has been constructed in [10] for infinite-dimensional 
representations of Verma module type. The spins s and i are assumed to be generic. The case of (half)- 
integer spins demands an additional refinement. Indeed, the limit s —^ in Eq. ([5|) has to be calculated 
carefully, since the divergences arise in both factors. In [8] it has been shown that at (half)-integer spin 
2s = re € Z>o the operator ([5|) can be restricted to a finite-dimensional invariant subspace in the first 
space of the tensor product. The restricted operator acts on a tensor product of the (re -|- l)-dimensional 
space (where spin s = ^ representation is realized) and an infinite-dimensional space (where spin £ 
representation is realized). In other words it is a (re -|- 1) x (re -|- 1) matrix, whose entries are differential 
operators acting on the space of polynomials C[z\. In [8] the restriction of the R-operator has been 
calculated and the generating formula for its matrix matrix elements has been found. More exactly, the 
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R-operator being applied to the generating function {zi — x)'^ of the finite-dimensional representation 
in the first space and to a polynomial ‘h(z) from the second spac^ gives the following 


- xY ^{z) = (6) 

= {z- {zi - z)“+t+^+i (zi - -^-1 (z - $(z). 

Expanding both sides of Eq. Q in powers of the auxiliary parameter x we recover matrix elements of 
IRi2(u|^,£). If we choose the second spin in Eq. ([6]) to be (half)-integer as well 2 £ = m G Z>o, then 
we immediately obtain the restriction of the R-operator to a (m -|- l)-dimensional representation in the 
second space. Indeed, applying Mi2(u|^, to the generating function <h(z) = {z — y)™ of the finite¬ 
dimensional representation in the second space, we find the generating function for matrix elements of 
1^12 y), which is a (n-|-l)(m-|-l) x (n-|-l)(m-|-l) matrix solving the Yang-Baxter equation, Eq. (fTj). 

The formula ([6|) contains in a compact form all matrix elements of the restricted R-operator. However 
the matrix form of the restricted R-operator is still rather obscure. Using the formula ([6]) as a starting 
point, we will infer an explicit formula for the restricted R-operator as a matrix of differential operators. 
Moreover we will see that this matrix is organized very simply and that it is much more transparent 
than In order to get accustomed to ([6]) let us consider several examples. 

In the case of restriction to two-dimensional space (spin s = ^) the formula ([6]) gives rise to the 
quantum L-operator [ 15 ]. In order to see it, let us choose the following basis in C^: ei = zi, 62 = 1 . 
In matrix notations ei = ( 1 , 0 ), 62 = ( 0 , 1 ). Equating coefficients by equal powers of the auxiliary 
parameter x in both sides of 

1Ki2('W — ^ 1 ^)-^) {zi — x) $(z) = {z — x)““+^+^ (zi — 2;)“+^+^ dz (zi — 2)“““^ {z — xY~^ d>(2;) 


yields the action of the R-operator on the basis elements 61,02 

l^i2(^i — ^1^)^) ei = 61 {zd — i + u) + e2 {—z^d + 2£z ), 

Ki2(w - ||^,^) 62 = 61 (9 -I- 62 (u -I- ^ - zd). 

We tacitly assume that both sides of the previous equalities are applied to an arbitrary polynomial <^(2:). 
Thus the matrix of the operator Mi2(u — |||)^) in the chosen basis is the following 


^12{u - l\lJ) 


I' u — £ + zd d + S S \ 

\—z‘^d + 2 £z u + £ — zd) V®"*" u — Sj' 


( 7 ) 


It does coincide with the L-operator. The implemented shift of the spectral parameter simplifies the 
previous formula. A straightforward calculation enables to check that the L-operator is a product of 
several upper-triangular and lower-triangular matrices 


Ki2('W - ^\IJ) 


( 1 0\ A 0 \ A (ui 0\ A 0\ 
V VO W VO V VO V V^ V 


( 8 ) 


where we introduce linear combinations of the spin and the spectral parameter 


ui = u — £ — 1 , U2 ^ u + £ . 

The factorization formula ([8|) is rather natural, since both the initial infinite-dimensional R-operator, 
Eq. dSj), and its restriction, Eq. Q, have the factorized form. This leads to a reasonable question: 

® In order to simplify notations we use z instead of 22 . 

Inessential normalization factors in and m are different. 
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does there exist a factorized matrix form of the R-operator at (half)-integer spin s = ^ which would be 
analogous to ([5]) ? 

In order to answer this question let us consider a more intricate example of spin s = 1, i.e. the 
restriction to (n = 2 ). Straightforward application of the formula ([U]) yields the following matrix of 
the operator Mi 2 (u — written in the basis ei = zf, e 2 = zj, = 1 

n2 + n(2S-l) + S(S-l) -nS" + SS" (S-)^ 

-2wS+ + 2S+S -u-2S'^+£{£ + !) -2nS-- 2S-S 

(S+)2 -US+-SS+ u2_^(2S + 1) + S(S + 1) 

The previous matrix entries are represented through the algebra generators, Eq. ([3j). It is easy to check 
that the matrix can be decomposed in a product of several more simply organized triangular matrices 

/I 0 0 \ /I 0 0 \ /I d d‘^\ /ui{ui-l) 0 0 \ /I 0 0 \ 

{-2z 1 0 0 U 2 - 1 0 0 1 2(9 0 m - 1 0 2z 1 0 . (9) 

\ -z 1/ \0 0 U2{U2-1)/ \0 0 I J \ 0 0 ij \z^ z 1/ 



Considering more examples we are able to guess the factorization formula for the matrix of the operator 

- , - V I Itl I tlU ItAJUIW <J ^ 1 

2 I 2 • 


{n-my 


written in the basis ei = z^ ,e 2 = z^ , • • • , e„_|_i = 1 

Ki2(w- §!§,£) = Z-^U+{u2)DU-{ui)Z . 

For the few first triangular matrices Z and D (at spin s = ^ , 1, |, 2, • • •) we obtain 


( 10 ) 


^(|) “ 


^ih) = 
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2 ^ 1 0 
’ 2 1 , 


•^(|) - 
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29 ,79(1) = 

VO 

0 

1 / 


/I 0 0 0\ 

3z 1 0 0 

3^2 2 z 1 0 


0 1 2 d 302 
0 0 1 3(9 

Vo 0 0 1 / 


, Z(2) = 


, D(2) = 


/ 1 


0 

0 0 

o\ 


4z 


1 

0 0 

0 


6^2 


3z 

1 0 

0 


Az^ 

3^2 

2 z 1 

0 


\z^ 


z^ 

2;2 z 

1/ 


n 

d 

92 

93 

94 

\ 

0 

1 

29 

392 

493 

0 

0 

1 

39 

692 

0 

0 

0 

1 

49 


VO 

0 

0 

0 

1 
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( 11 ) 


( 12 ) 

Thus one infers immediately the general pattern. The diagonal matrices U~^{u) for the same values of 
the spin are the following 


C/^V) = diag(l,u) , = diag(l,'u - l,u('u - 1 )), 

= diag(l,u -2,(u- l)(u - 2),u{u - l)(u - 2 )), 

^( 2 ) ~ diag(l, u — 3, (u — 2){u — 3), {u — l)(?x — 2){u — 3),u{u — l){u — 2)(n — 3)), ■ ■ ■ (13) 

The eigenvalues of U~{u) are in a reverse order 


17(1) = diag('u, 1 ) , C/(^) = diag(u(M -l),u- 1 , 1 ), 

C/(7) = diag('u(M - l)(ri - 2), {u - l){u -2 ),u- 2,1), 

17(2) = diag(ri(u — l)(ri — 2){u — 3), (u — l){u — 2){u — 3), {u — 2)(rt — 3), u — 3,1) • • • (14) 
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The examined examples lead to a transparent factorization pattern expressed by Eq. (llOp . The fac¬ 
torization formula m offers considerably more explicit description of finite-dimensional solutions of 
the Yang-Baxter equation, Eq. m, as compared with all other known approaches. The factorization 
formula (I1U|) is equivalent to the generating formula for matrix elements of the restricted R-operator, 
Eq. ([6]), hence it confirms the efficiency of the result ([6]) obtained initially in [8]. Let us mention that 
it would be difficult to guess the factorization formula (jlOh taking as a starting point the classical 
result dH). 

The quantum Lax operator (also called L-operator) is a 2 x 2 matrix (the fundamental representation 
of a rank 1 algebra is two-dimensional), and its matrix entries are linear in generators of the symmetry 
algebra. The L-operator can be factorized in a product of several more simple matrices in the case of 
si 2 symmetry algebra as well as in the case of its trigonometric and elliptic deformations [Hill 1(125] . This 
observation helps a lot in solving the RLL-relation, which imposes severe constraints on the infinite¬ 
dimensional R-operator [Tl llllfTH] and eventually enables to find the general solution of the Yang-Baxter 
equation, Eq. ©• The purpose of this note is to show that more general solutions of Eq. m than the 
L-operator can be factorized as well. 

In the next Sect, we will prove the factorization formula (jlOP - Besides the reduction formula ([6]) 
for si 2 the analogous restrictions (to a finite-dimensional representation) of the general R-operator has 
been obtained in the paper [8] for the Lie group SL(2,C) and for the modular double of Uq{s(. 2 ) [TT] - 
In the accompanying paper [9] the analogous restriction of the general R-operator m has been carried 
out for elliptic deformations of si 2 , which are the Sklyanin algebra and the elliptic modular double [38] . 
The matrix factorization of SL(2, C)-symmetric R-operators does not differ essentially from the formula 
m, since finite-dimensional representations of SL(2, C) are tensor products of two si 2 representations. 
So we will not consider factorization for SL(2,C). In Sect. [3] and U] we deal with the modular double 
and obtain counterparts of the results presented in the current Sect. The trigonometric factorization is 
given by the formula (1351) . In Sect. [5] we consider elliptic deformations and find the elliptic factorization 
formula (IHHp . 


2 Rational factorization 


In this Sect, we will prove the matrix factorization formula for the restricted s£2-iiivariant R-operator, 
Eq. (jinj) . The reduction formula ([6]) obtained in [8] will be a starting point for us. The proof consists 
of two steps. Firstly, we rewrite the matrix formula (lIOp in an operator form. Secondly, we act by this 
operator on the polynomial (zi — x)^ ^{z) and transform the result to the form ([6]). 

For the sake of simplicity let us consider firstly an example such that the spin s = 1. We are going 
to rewrite the factorization formula (lin|) at s = 1 in an operator form. Recall that the matrix of the 
operator Mi2(ri — 1|1,-^) does factorize, Eq. Q. It is constructed out of matrices (see 

Eqs. (HH), (HU), (HHI), (HID) 

Mi2(n - f |I,.^) = C/(+)(n2)L>(i) t/(7)(ni) Z(i) . (15) 

Each matrix factor in the previous formula has an operator counterpart. The matrices and 
have a simple exponential form 


Z(i) = exp (2 ;D(i)) ; D(i) = C(i) exp (5D(i)) C(i) , 

where we introduce the numerical matrices 
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0\ 
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A 
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; C(i) = 0 

1 

0 . 

(16) 

\0 

1 

0/ 

Vi 

0 

0/ 
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The matrices and are related to each other by the similarity transformation 


hi) 


[/+(n2) = C(i)C/(-)(«2)C(i). 


(1)^ 


Substituting the previous expressions in Eq. m and taking into account that C(i) = 11 we rewrite 
the matrix (jl5l) as follows 

Mi 2 (u - = exp (- 2 ;D(i)) C(i) (^ 2 ) exp («9D(i)) C(i) UZJui) exp (^D(i)) . 


The matrices that constitute the previous expression are matrices of some operators in the basis ei = zf, 
e2 = zi , 63 = 1: 

• the lower-triangular matrix D(i), Eq. (|16p . is a matrix of the differential operator dz^ with respect 
to the given basis; 

• the matrix Eq. (|16l) . is a matrix of the inversion operator Ci = C (8) 11 : —>■ Zi~^ at n = 2 

with respect to the given basis; 

• the diagonal matrix (u) is a matrix of the operator at n = 2 with respect to the 

given basis. 

Straightforward generalization of the considered example s = 1 (i.e. re = 2) to arbitrary re yields an 
equivalent operator form (with respect to the basis ei = z” , 02 = z^~^ , • • • , e^+i = 1) of the matrix 
factorization formula, Eq. m, 

« / n\n r o \ CS ^{zidi + U 2 + l-n) r(2;i5i+rei-k 1 - re) 

Mi 2 (re-f| 2 ,£) =exp(-z5i) Ci- + 1 - n) - exp(59i) Ci- + 1 - n) - 

(17) 

We proceed to the second step of the proof. We are going to act by the operator (flTll on the polynomial 
{zi — x)”<f>(z) and to verify that the result does coincide with the formula for matrix elements of the 
R-operator, Eq. ([6|). Thus we act sequentially on {zi — x)” <h(z) by the operator factors from Eq. (fT7|) . 
At the first step we perform the shift zi —)■ zi -h z, 

n I 

exp (z di) (zi - x)*^ $( 2 ) = (zi + z - xY d>(z) = ^ 


fe=o 


At the second step we act by the operator on the previous expression. The action of this 

operator on z\ is equivalent to the substitution z\d\ ^ k, 


E 

k=0 


re! r(A:-I-re 2 + 1 — re) ^ 


k\{n — k)\ r(re 2 + 1 — re) 


zf (z - $(z). 


(19) 


At the third step we act by the operators Ci and exp (ddi) on z^ that is present in Eq. (flUD . 

( — til 

exp (aai) Cl 4 = exp (5^1) zr^ = (zi + 5)"-^= = u _ k _ u ^ ’ (^0) 


p=o P^-{n-k-p)\ 

and at the last step we apply exp {—zdi) Ci that is present in Eq. ([20]) 

'{U2 + 

r(re2 


exp (-. a) c, = (=. - 


r(re 2 + 1 - re) 
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Finally, gathering the previous expressions, Eqs. (fT9l) . (pOl) . (l2T]l . we obtain 


Mi2(n-§!§,£) izi-xr^z) = 


E 

k=0 


ni 


k\{n — k)\ r(tti 


r(m + l-n + fc) ^ {n-k)\ Tiu2 + l-k-p) _ .k+p _ .n-k . 
r(ui + l-n) ^^pl{n-k-p)] T{u2 + l-n) 


— n) 


( 22 ) 

The previous formula is an operator reformulation of the matrix formula (|10l) . More exactly, according 
to Eq. (I22p . the matrix (llOh is applied to the (n + l)-dimensional vector {zi — x)"' and operator entries 
of the matrix (jlOp act on a polynomial In order to complete the proof we need to show that the 

right hand side of Eq. ()22p coincides with the right hand side of Eq. ([6]) where the spectral parameter 
is shifted u 


u — 


2 ’ 


Ki2(w - {Zi - xY ^{z) = (z- (zi - (zi - z)-“2-i+n (2 _ _ ^ 23 ) 

The proof of the needed identity will be based on the Cauchy’s differentiation formula 




dA 


(A - z)P+^ 


F(A) 


(24) 


where the closed contour around A = z does not encircle singularities of an analytic function F(A). 

Let us consider the right hand side of Eq. (1221) . We reduce fractions canceling {n — k)\ and change 
the summation index in the second sum p ^ n — k — p 


n—k 


n! r(ni + 1 — n + A:) 

^ Vin-i -1-1 — rj'\ 


1 


A! r(ui + 1 — n) p!(^ “ ^ “ P)' r(tt 2 + l —n) 


T{u 2 + 1 n + p) ^ 


Eurther we exploit the integral representation 

n—k 


for the underlined factor. 


n! r(rii + 1 — n + k) 

^ 'Pfni -1-1 — n'\ 


r(ui + l —n) ^^p\{n — k — p)\ T{u2 + l — n) 


1 r{u2 + l-n + p) (n-fc-p)! ^^ 


27ri 


dAiA-^^ 

(A - 5 ;)n-fc-p+i 1 1 27ri 


dA 


-{zi-zY{A-xYx 


E 


r(ni + 1 — n + A) f A — z 


fc_o ^!r(ni + l-n) \A-xJ ^ p!r(tt2 + l-n) 


k n-k 

E 


(A - zY+^ 

r(tt 2 + l — n + p) / A — z^^ 


Zl — z 


^A) 


(25) 


Both summations in the previous formula can be extended freely to infinity. Indeed, the unwanted 
terms contain (A —z)”*, m > n + 1, and, consequently, disappear being integrated over A. The emerging 
power series are of binomial type 




k=0 


k\T{a) 


so the sums in Eq. 


can be evaluated explicitly 


n\ 


dA 


27ri 7 {A-zY+^ 


{zi - zY (A - xY 


1 - 


A — z 
A — X 


n—ui — 1 


1 - 


A — z 
Zl — z 


n—U 2 —l 


cb(A) 
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Further we rearrange the factors in the previous expression and calculate the contour integral according 
to Eq. (0^ that produces immediately the desired result (E5|) . 


(z - (z, - (^1 - = 

= {z- {zi - + l {zi - z)-“2-l+n _ 3.)«l + l _ 

Thus the identity (12311 along with the matrix factorization formula for the operator Mi 2 (u — 
Eq. (fTOl) . are proven. 


3 Modular double 


In this Sect, we consider solutions of the Yang-Baxter equation, Eq. that are invariant with respect 
to the modular double. The modular double of the quantum algebra Uq{s(. 2 ) has been introduced by 
Ludwig Faddeev in m- This algebra is formed by two sets of generators E , K , F and E , F , K. The 
standard commutation relations for the generators E,K,F, which form the quantum algebra Uq{si 2 ) 
with the deformation parameter q = (we assume that r G C\Q, i.e. q is not a root of unity) 

[E, F] = , KE = gEK, KF = q-^FK, (26) 

are supplemented by analogous commutation relations for E, F, K with the deformation parameter 
q = In addition, the generators E and F commute with E and F; the generator K anticommutes 

with E and F; K anticommutes with E and F; K and K commute. 

The representation theory of the modular double has been elaborated in a number of papers, see 
for example [UIITIIIHIEIIEI] and references therein. We will use the following parametrization t = ^, 
where w, a;' G C, Imw > 0, Imw' > 0, are constrained by ujuj' = —j- Then 

q = exp [muj'/uj) , q = exp [iTruj/uj') , 

so the change q ^ q is equivalent to u ^ uj'. We will also profit from the notation lu" = u; + ui'. 

Further we deal with realization of the modular double generators by finite-difference operators 
Ks = 7rs(K),Es = 7rs(E),Fs = 7rs(F) acting on the space of entire functions rapidly decaying at 
infinity along contours parallel to the real axis. The representation vr^ is parametrized by a complex 
number s, which we call spin. The generators have the following explicit form 


K, = e-^P 


(g-r')E. 

(g-r')F. 


ITTX 

= e “ 




_ ITTX 

= e “ 


tiP+s+^'') 


e2ui 


- e-^^P+^+^'') 


(27) 


where p is a momentum operator in the coordinate representation, p = The formulae for 

generators , E^ , F^ are obtained by a change cj ^ w' in Eq. ([27)1 . 


The noncompact quantum dilogarithm naturally arises in the representation theory of the modular 
double. In the context of quantum integrable systems it has been found first in |16] . The properties 
of this special function have been thoroughly examined in [181143] . We will need not the quantum 
dilogarithm itself but another closely related special function dehned by the integral 

+00 \ 

1 f dt sm{at) cos{zt) \ 

2 J t sin(a;t) sm{oj't) j ’ 

— OO / 
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Da{z) = exp 


(28) 









where the contour goes above the singularity at t = 0. The R-matrix of the Faddeev-Volkov model is 
expressed through this function [51[5D]. A number of identities for the T)-function are contained in [B], 
It naturally arises as an intertwining operator of equivalent representations of the modular double. 

Further we indicate some basic properties of the D-function that we will need. The function Da{z) 
is even, 

Da{z) = Da{-z) ; Da{z)D-a{z) = 1 ; Do{z) = 1. (29) 

From the definition (1281) we infer that it is symmetric with respect to the change uo oo'. The D-function 

satisfies a pair of finite-difference equations of the first order 

Da{z-u:') ^ cos^{z - a) _ Da{z - u) ^ cos^{z-a) 

Da{z + uj') cos^{z + a) ’ Da{z + u) cos^{z + a) 

Consequently 2a; and 2uj' have the meaning of its quasi-periods. 

At generic spin s the representation vr* is irreducible and infinite-dimensional. However it is not of 
a Verma module type, since the representation space does not contain a highest-weight vector n(a;); 
Fsn(x) = 0, F 5 n(x) = 0, Ksfl(x) = An(x), Ksn(x) = An(x). The representation is a deformed 
analogue of the principal series representations of the Lie group SL(2,C). The situation drastically 
changes at spin s = Sn,m. = — w" — nu — muj', where integers n,m € Z>o enumerate points of a quarter- 
infinite lattice on the complex plane (or a line, for real u/oj'). In this case the representation is 

reducible, and a (n -|- l)(m -|- l)-dimensional irreducible representation decouples. Since we will need 
finite-dimensional representations, let us consider their structure in more details. The basis of the 
(n -|- I)(m -|- l)-dimensional representation is formed by monomials 

xn-2kj^m-2i A: = 0,1, • • • , n , / = 0,l,---,m, 

with respect to the variables 

X = X(x) = , X = X(x) = e^^. (31) 

Therefore any finite-dimensional representation of the modular double is a tensor product of finite¬ 
dimensional representations of Uq{si 2 ) and Uq{si 2 )- For our purposes finite-dimensional representations 
of spin s = Sm = —uj" — mu', m G Z>o, will be sufficient. Thus we will deal with only a half of the 
modular double. By means of the D-function, Eq. (1281) . the basis vectors I = 0,1, ■ ■ ■ ,m, can 

be arranged in a sole function. Indeed, Dmu)>{x — y) is a generating function of the (m -|- l)-dimensional 
representation. It reduces to a linear combination of exponents by means of certain contiguous relations 
similar to (|30|) . 

m—1 

Dmco'{x-y)=ll(Y-^Xq^-^ + YX-\-^+^) , Y = Y{y) = e^y, (32) 

1=0 

where y is an auxiliary parameter. 

4 Trigonometric factorization 

After a brief introduction to the modular double given in the previous Sect, we will consider solutions 
of the Yang-Baxter equation, Eq. that are invariant with respect to this quantum algebra. The 
invariance means that R-operators commute with the co-product of all six generators. Let us note 
that the symmetry restriction imposed by Uq{si 2 ) is not sufficient to fix completely (up to an inessential 
normalization factor) the R-operator. Initially the R-operator for the modular double acting on a tensor 
product of two infinite-dimensional representations (8> was constructed in [6] in a form similar 
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to Eq. Q where the role of the Euler beta function is played by the -D-function, Eq. (1281) . The remarks 
made in Sect. [T] concerning the formula Q are equally valid for this realization of the R-operator for the 
modular double. In [7] relying on the realization (|27p of the algebra generators a more explicit formula 
for the R-operator has been proposed. There we demonstrated that the R-operator can be realized 
as an integral operator, and it factorizes in a product of four Eaddeev-Volkov type R-matrices [20] . 
In [29[l3nj an explicit hyper geometric formula for the R-matrix of Uq{si 2 ) acting on a tensor product of 
two highest-weight representations has been discovered. In [24] a universal factorization formula for the 
trigonometric R-operator has been obtained by means of the universal R-matrix. 

In [ 8 ] the integral R-operator for the modular double acting on a tensor product of two infinite¬ 
dimensional representations ^ (8* vr^ has been used as a tool to produce finite-dimensional solutions 
of the Yang-Baxter equation, Eq. ©• There the restriction of the R-operator to a finite-dimensional 
representation in the first space at spin Sn,m = —w" — nu — mu', n, m ^ Z>q , has been found. Let us 
mention that in m a similar restriction has been implemented for Racah-Wigner 6 j-symbols. 

As we have already explained in the precious Sect., we are going deal with finite-dimensional repre¬ 
sentations only of spin Sm = —mu' — u". A generic finite-dimensional representation of spin Sm,n can 
be taken into account as well without considerable changes in the following reasoning. The R-operator 
acts on the generating function of finite-dimensional representations at spin Sm, Eq. (I32p . according to 
the following formulqj 


Kl2(w|sm,s) • Drai^i{xi^)^{x2) = Du^{xi2)^ 

{^ 23 ') ■ ^muj' {P 2 ') ■ {^ 12 ) 7)iii(3J23) d?(x2) , 


(33) 


where X 3 is an auxiliary parameter of the generating function. We recall the shorthand notation Xij = 
Xi — Xj. Instead of the spectral parameter and spin s we prefer another pair of parameters 

s s 

ui=u + U2 = u- 

In these variables the final result will take a more simple form. Dmui'{p 2 ) from Eq. (|33p is a finite- 
difference operator. It factorizes in a product of m finite-difference operators of the first order due to 
the formula (1321) . 

The generating formula (l33]) uniquely specifies the solution of the Yang-Baxter equation, Eq. (flj). 
that is a (m -|- 1) x (m -|- 1) matrix with operator entries. According to Eq. (|33p . the entries are finite- 
difference operators of order m. The formula (|33l) enables to obtain a more explicit realization of the 
restricted R-operator. With respect to the basis 

j = l,...,m + l, (34) 

where Xi = Yi(xi), Eq. (|3ip . the matrix factorization formula for the restricted R-operator is valid 

1 ^ 12 ('u|Sm, s) = Z M{u 2 ) D M (tti) Z~^. (35) 


The previous formula is a trigonometric counterpart of the rational factorization, Eq. ()10ll . The matrices 
Z and D are diagonal 


{Z)kj = 6kj vf; {D)kj = 6kj q(m-l){m+2-2k) ^(m+ 2 - 2 k)u ^'82 ^ 


(36) 


The coordinate X 2 is present only in the matrix Z, and the momentum operator p 2 is present only in 
the matrix D. The numerical matrix M{u) is given by the following hyper geometric sum 




( 9 ^;9^)fc-p-i(e^; 9^}m+2-j-p+p 


jj2{2p-j-k+2)+m ^ 


^Here we implemented a shift of the spectral parameter u as compared with [8]. 
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where the summation over integer p is from max(0, k-\-j — 2 — m) to min(A: — 1, j — 1); the q-Pochhammer 
symbol {q^]q^)k = (1 - - q^) ■ ■ ■ {I - q^^)] U = U{u) = e^. 

In order to illustrate the formula (I37h we indicate the first few matrices M{u), m = 1,2, 3. Here we 
use shorthand notations + m), 





m(2) 


/ [/2 1 C/-2 \ 

[q + q-^ gC/2 + g-lC/-2 g + g-l 1 

\ 17-2 ^ ^2 j 


/ t/3 U U-^ 17-3 \ 

. ,( 3 ) _ {q^ + l + q~‘^)U q^U'^ + {I + q~‘^)U~^ q~‘^U~‘^ + {I + q^)U {q^ + I + q~‘^)U~^ 

(g2 + 1 + g-2)c/-i + + + {I + q-‘^)U-^ {q^ + l + q-‘^)U ' 

\ [/-3 17-3 ^ J 

In the case m = 1 the R-operator being restricted to the fundamental representation turns into the 
quantum Lax operator [6]. The factorization of the L-operator of the XXZ spin chain has been discovered 
firstly in [3] in the context of the chiral Potts models. 

The rest of this Sect, is devoted to the proof of the trigonometric factorization formula, Eq. (I35|) . 
We rewrite the finite-difference operator D^i{p 2 ) from Eq. ([33]l as a sum of shift operators 


m+1 

Dm.'{p 2 ) = Yl dj ^ ^ 3 g) 

i=i 

An explicit expression for numerical coefficients dj will not be relevant for a while (see Eq. (I47h ). Then 
we rearrange the factors in Eq. (|33p . We collect all functions depending on U 2 to the left of the shift 
operators and all functions depending on ui to the right of the shift operators, 

m+1 

Kl2(w|Sm,s) • Dm.'{xi'i)^{x2) = E dj Du2{^12) ii2+maj^ {xi 2 + {2j - m- 2)ujjx 

i=i 

• e(™+2-2jV'(92 Du,, {X 23 ) D_u,^+m.'{x23 + (2j - 171 - 2)w') ^{X 2 ) ■ (39) 


So the coordinates are present in Eq. (f3^ only in the form of the function Du{x)D-u+mui'ix + {2j — 
m — 2)oj'), where j = l,...,m + l. By means of contiguous relations similar to (|30p one can check that 
this function is given by the following finite product 


Dui^x')!)—u-\-mw'ix ~\~ {2j TTl 2)(U ) — 

m—j j—2 

= 17”^+2-2i gi-f-1 X™ (l -7 gX-317-3 ^2fc^ JJ (^ + X-^Uq^^^ . (40) 

k=Q k=0 

Expanding the right hand side of Eq. (IdOp we obtain the following sum 

m+1 

Duix)D_u+m.'ix + ( 2 j -m- 2)uj') = ^ djk{u) , (41) 

k=l 

where djk{u) are some numerical coefficients, which will be calculated afterwards (see Eq. (l48p L 

Now we are ready to calculate the matrix of the operator Mi 2 (M|sm, s) with respect to the basis (fMp . 
We substitute the expansion (l3Tp into Eq. ([5911 to the left and to the right of the shift operators. 
Then we expand both sides of Eq. (I3^ in powers of X 3 = ^ 3 (^ 3 ), Eq. ([3T]) . The generating function 
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can be expanded with respect to the basis (IMt) and simultaneously in powers of according 
to Eq. (1521) . Equating coefficients by powers of in both sides of Eq. (I55|) yields 


^I2(w|sm,s) • dkek^{x2) = 

/ m+1 


m+l 

i=l 


L 


dji{u2) dj e("^+2-2i)^'92 


djk{ui)X^ 


m-\-2—2k 


$(rE2). 


(42) 


Matrix entries of the operator M are coefficients in the expansion of the vector Me^ with respect to 
the basis ([MI): 1^6^ = (M)ik- Consequently the formula (1421) produces immediately the matrix 

entries of Mi 2 (u|sm, s)j 

(Mi 2 (n|.^, s)),k = dp,{u2) dp) ^^+2-2k'^ ( 43 ) 

In the previous formula we tacitly imply summation over repeated indices p, j; the matrix entries are 
presented in the operator form, and we omit an arbitrary function <h(x 2 ). Thus the right hand side of 
Eq. (1431) is factorized. It is a product of three matrices: the diagonal matrix containing shift operators is 
sandwiched between two matrices made of numerical coefficients djk{u), dk and coordinates . 

Stripping off the diagonal matrices that contain X 2 from the lateral matrices we obtain the factorization 
formula 

Kl 2 («|'Sm, s) = Z M 2 {u 2 ) D Mi{ui) Z~^ . (44) 

The diagonal matrix Z is defined in Eq. (j36p . The diagonal matrix D is slightly different from D defined 
in Eq. ([36]) . 

The numerical matrices Mi{u), M 2 {u) are constructed out of the expansion coefficients djk{u), d^, 

(Mi(u)),, = ^^ ; {M2{u)\, = dkdki{u). (45) 

The factorization formula ()44p is slightly different from Eq. (j35p . In order to recast the formula (I44p 
into (1351) . first of all we need to find the coefficients djk{u), dk defined by expansions ([381) and ([Ml) . This 
goal is easily accomplished by means of the q-binomial theorem 


m—1 




A:=0 


xq"^^ \ = 


E 

A:=0 


id 




(g2;g2 


j(q' 2 , q‘^')Yn—k 


(46) 


Indeed, the function Dmu}/, which produces coefficients dj, Eq. (l38p . is just the product ([32]) of the 
type (]46p . Consequently, 

, (g 2 ;g 2 )^^(fc-l)(fc-m-l) 

iq^;q^)k-i{q^;q^)m-k+l' ^ ^ 

The coefficients djk{u), Eq. (|lTp . can be extracted from the product of two q-binomial sums. We omit 
the details of the calculation that results in 


djk{'^) 


E 


( 9 ^; q^)k-p-i{q^] q^)m-j+2-k+p 


jj2{2p—j—k+2)+m 


(48) 


The summation limits over integer p in the previous formula are the same as in Eq. ()37p . 
Let us define 

djk{u) = dj djk{u) 2 . 
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Substituting the explicit expressions for the coefficients (H71) . (H8|) in the definition of djk{u), one straight¬ 
forwardly checks that it is symmetric in indices j,k: djk{u) = dkj{u). This observation enables us to 
simplify Eq. (HU) . We separate the diagonal matrix 6ik d ^, move it from M 2 (u) towards Mi (u), Eq. ()35]) , 
and cancel it. Thus instead of a pair of different matrices Mi(u) and M 2 {u) a pair of identical matrices 

1 A I 1 7Tl(7Tl+l) _ _ _ _ . „ „ . _ _ _ 

{M{u))j^j = djk{u) 2 ^ Eq. ([37]), is present in the factorization formula (13511 . The formula 

(f35]l is proven. 

It would be interesting to relate the factorization formula ()35|) to explicit expressions for R-matrices 
from [291130] as well as to the universal factorization formula |24j . 


5 Sklyanin algebra and elliptic factorization 

In this Sect, we will factorize solutions of the Yang-Baxter equation, Eq. whose symmetry is 
encoded by the Sklyanin algebra [36]. The Sklyain algebra is a two-parametric deformation of si 2 or 
a one-parametric deformation of Uq{si 2 )- It serves as a dynamical symmetry algebra of the 8 -vertex 
model [T]. The four generators S^, S^, of the algebra respect commutation relations 

S" S“ = i • (S° S°) , 

gO s“ - S“ S° = i ■ (S^ S^) , (49) 

where the triple {a,l3,'y) is an arbitrary cyclic permutation of (1,2,3). The structure constants Ja/s = 
'j ^ a, P, are expressed through the Jacobi theta functions (we assume t/ € C and 6 a{r]) 7 ^ 0, a = 
1,...,4) 

Ji = 02 ( 2 ??) 02 (O) 02 “'(^) ; J 2 = es{ 2 rj) 6 s{ 0 ) 6 f{rj) ; J 3 = 6 ^{ 2 rj) 6 ^{ 0 ) 6 f{rj). (50) 

We adopt shorthand notations 6 a{z\T) = 6 a{z), a = 1, ■■ ■ , 4, for theta functions with modular parameter 
r G C, Im(r) > 0, 

0 i(z|r) = ei{z) = -Y, . 

The rest three theta-functions are obtained by shifts of the argument of 61 by quasi-period halves 
92{z\t) = Oiiz + 1 \t) , 93{z\t) = e~~^''^^92{z + ^\t) , 9i{z\T) = 93{z + \\t) . 

Besides theta functions 9a{z) with modular parameter r we will need as well theta functions with 
modular parameter which we denote as follows 

h{z) = 93 {z\^) , 94^{z) = 94^{z\^) . (51) 

Two types of theta functions are related to each other by the identity 

2 9i{x + y)9i{x -y) = 9i{x)h{y) - 9A{y)h{x) ■ (52) 


In the previous Sect, we have seen that the noncompact quantum dilogarithm (more exactly, the 
D-function) is omnipresent when one deals with representations of the modular double. In the case of 
elliptic deformation the same role is played by the elliptic gamma function [35] 


T{z) = T{z\T,2y)= n 


n,m=0 


_ g—27ri2^n+l^m+l 
_ Q27viZpnqm 


p = e 


27riT 


q = 


(53) 
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where |p|, 1^1 < 1- This function possess a number of remarkable properties. We will need the reflection 
formula 

r(z) r(— z + 2?7 + r) = 1 (54) 

and its quasi-periodicity at the shift by 2??, 

_ i_ 

T{z + 2rj)=R{T)e^^^9i{z\T)r{z) , R{t) =. (55) 

^\Pi Pjoo 


In the following we extensively use the short-hand notation r(± 2 ; ± x) := r(z -|- x)T{z — x)T{—z -|- 
x)T[—z — x). Various connections between the Sklyanin algebra and elliptic hypergeometric functions 
were considered in [32H341I38] . 

Let us briefly outline some basic facts about representations of the Sklyanin algebra. It admits a 
highly nontrivial explicit realization of generators as first order finite-difference operators with elliptic 
coefficients found by Sklyanin in his pioneering paper 133, 

ei{2z) 


9a+i {2z- g + rj) - Oa+i {-2z - g + rj) e 


Vv 


(56) 


The operators depend on a parameter g (z C called the spin. They act on the space of holomorphic 
functions of z. In Eq. (j56p we use unconventional similarity transformation by means of whose 

meaning is explained in m- At generic g the representation ()56l) is infinite-dimensional and irreducible. 
However, for a discrete set of spin values g = g^ ^ {n + l)r] n G Z>o, a (n -|- l)-dimensional 
representation decouples. The finite-dimensional representation can be realized in the space 0^ of even 
theta functions of order 2n. It is formed by holomorphic functions that are even f{z) = f{—z) and have 
simple quasi-periodicity properties under the shifts of z by 1 and r: 


f{z + 1) = fiz) , f{z + t)= e-2---4--/(^). 


The action of the generators (j56h at spin g = g^ is invariant and irreducible on this space. One can 
easily check that the monomials constructed out of theta functions (fCT) form a basis in 

the space 

= [h{z)Y [ 04 :{z)Y~^ , j = 0,l,---,n. (57) 

The elliptic gamma function, Eq. (|53l) . enables to combine the basis elements into a sole 

object. Indeed, T {^z x + gn) is a generating function of the (n -|- l)-dimensional representation, and 
it depends on an auxiliary parameter x. Owing to Eqs. (1521) . (fSTll . (1551) it reduces to a product of linear 
combinations of O^iz) and 04 (z), 


n—1 

c-T {^zz^x + gn) = [h{z)h{x + {n-l- 2r)r]) + O^iz) 03 (x -h (n - 1 - 2r)g) ] , (58) 

r=0 


where an inessential numerical constant c = (—2)"'R“^"'(r)e“~"'. The previous product is equivalent 
to a linear combination of basis vectors Eq. (j57p . with some coefficients 'tjj^'^\x) depending on 

the auxiliary parameter x. The generating function, Eq. (I58p . is invariant under the change z ^ x, 
hence it contains the second natural basis (^)}j=i, 


n+1 n+l 

T{^ZZfx + gn) = Y^ (^) = ^ PnU-jix) V’, 

i=i i=i 


(n), 
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(59) 


that is formed by products of O^lz) and 04 (z) with shifted arguments, 

n—1 

= Sym 0a, (2 + (n - 1 - 2r)r/) , o, € {3,4} , j = 0,1, • • • , n 

r=0 

where 03 appears j times and 04 appears n — j times; symmetrization Sym is over indices {a,}. 

Let us denote a pair of basis (1571) . (I5^ of the (n +l)-dimensional space 0^^ by and , 

ej = , fj = , j = 1,2, • • • ,n + 1. 

At n = 1 the representation is 2-dimensional, the bases coincide 

ei = fi = 04 ( 2 ;) , e 2 = f 2 = 03(^) ■ (60) 

At higher spins the bases are different. At n = 2 the representation is 3-dimensional, a pair of bases is 

ei = 0|(2;) , 62 = 04(^)03(5;) , e3 = O^iz) ; 

fi = 04 ( 2 ; - r])e 4 {z + r]) , f 2 = 04 ( 2 ; - r]) 63 iz + ??) + 03 (^ - ?/) 04 (^ + r/) , £3 = 03 ( 2 ; - r]) 93 {z + f]). 

These basic facts about the Sklyanin algebra and its representations will be sufficient for our purposes. 
Let us take a look at the corresponding solutions of the Yang-Baxter equation, Eq. &■ The symmetry 
restrictions imposed by the Sklyanin algebra do not allows to fix uniquely the solution Mi 2 (rt) acting on a 
tensor product of two infinite-dimensional representations specified by spins and < 7 ( 2 ) [IS]. However, 
more severe restrictions produced by the elliptic double enable to fix the R-operator unambiguously (up 
to an inessential constant). This R-operator has been constructed in |13] in a form of an integral operator 
acting on a tensor product of two arbitrary infinite-dimensional representations of the elliptic double. 
The integral kernel of this operator is given by the product of elliptic gamma functions, Eq. (|53l) . 
The proof that this integral operator with an elliptic hyper geometric kernel solves the Yang-Baxter 
equation is based on a number of sophisticated identities: the elliptic beta integral evolution formula |39l 
E], an integral Bailey lemma [40] . and elliptic Fourier transformation m- The elliptic beta integral 
evolution formula is equivalent to the star-triangle relation [3|. The elliptic double consists of two 
Sklyanin algebras, whose structure constants, Eq. (I50l) . are parametrized by 2r],T and T,2r], so their 
generators commute or anticommute with each other. Finite-dimensional representations of the modular 
double are equivalent (up to a sign) to a tensor product of finite-dimensional representations of the 
Sklyanin algebras. Since we are aimed at finite-dimensional representations and matrix realizations of 
R-operators, we will consider only one of two Sklyanin algebras constituting the elliptic double. 

In [9] the integral R-operator for the elliptic double has been taken as a starting point and restric¬ 
tions of this operator to finite-dimensional representations have been implemented. In particular, the 
restriction to a (n -|- l)-dimensional representation in the first space at spin < 7 ^ = (n -t- 1)77 + n £ Z>o, 
has been considered. The action of the R-operator on the generating function of the hnite-dimensional 
representation, Eq. (l58]l . is given by the formulel^ 


h 2 {u\gn , g) r(=F^l =F ^^3 + gn) ^(^) = 
r(=F^ TZ3- 2 


U _|_ £"±9' 


r(=F2i q: 2 ; - 1 - + 0 + 5) 


M(nr7) 


r(=F^l q: Z - I g^) 
r(=F^ T 23 - f + + 0 + i) 


4 >( 2 ). ( 61 ) 


An arbitrary holomorphic functions $( 2 ;) belongs to the second space where a generic spin g represen¬ 
tation is realized; 2:3 is an auxiliary parameter of the generating function. The finite-difference operator 
from the previous formula 

n 

Minri) = (z) ( 62 ) 

1=0 
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In order to simplify the notation we denote the complex variable zi by 2 . 
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is an intertwining operator of equivalent representations of the Sklyanin algebra. For the first time it 
has been constructed by A. Zabrodin in [33] . In m the factorized representation for the intertwiner 
has been found 

M{nrj) = Aa{nr] - rj) ■ ■ ■ Aa{r])Aa{0) ■ da"-(z) , o = 3,4. 

It is a product of finite-difference operators of the first order 


K{g) = e" " 


6i{2z\t) 


Oa{z + g + v) -Oaiz-g-g) e 


-Vdz 


e ^ 


The coefficients /3l'^\z), Eq. ([62]) . can be found in [T31I33] . 

Expanding (IbTTl by means of Eq. (158]) and equating coefficients on both sides of the formula that 
accompany the linear independent functions {<^^”^( 2 : 3 )} of the auxiliary parameter X 3 , we obtain the 
matrix form of the restricted R-operator 


Kl2(wbn , g)i^f'\zi) = ipf"\zi){Ri2{u\gn , 


with respect to the pair of bases (f57|l . {ej}^^l and , 


f^j = g^f'Hzi) , {j = il;f^\zi) , j = ,n-hl. 

The matrix elements are finite-difference operators of the n-th order whose coefficients are constructed 
out of theta functions. It happens that the matrix form of the restricted R-operator is more illustrative 
than Eq. dSIl). Indeed, this matrix solution of the Yang-Baxter equation can be factorized as follows 


, g) = Viui,z) D{z, d) C V^{U 2 , z) C . 


(63) 


The matrix D{z,d) is diagonal, and it is formed by the terms of the intertwining operator M.{nT]), 
Eq. (I62|), 

In the numerical matrix C only the antidiagonal is nonzero: (C);^ = 6 n+ 2 -i,j ■ We see that it is 
convenient to arrange the spectral parameter u and the spin g in the linear combinations 


ui 


u + g 
2 


, U2 = 


u- g 
2 


The matrix V consists of theta functions {V{u,z))ji = vj^\u, z) that are specified by the following 
defining relation 


n+1 n—l I 

vjp\z,u) = + z - u + ^ + 27?(| - / - r)) 01 (±x + z + u - ^ + 2g{^ - I + r)) . 

j=l r=0 r=2 

In view of Eq. (|52p . the function is a linear combination of theta functions 63 and 04 , whose 
arguments are shifted in a certain way. Each monomial contains j times 04 and n — j times 03 , i.e. 

Vjl"\z,u) = (-1)’*+!-^ Symj|^^^0„^_^ [±x + z + u - -b2r/(§ - l + r))x 

• n”=o + ^ - ^ + Y + 20(§ -l-r)), 

where G {3,4}. Let us note that an immediate corollary of the defining relation is vjj^\—z,u) = 
Vj^_^ 2 -iiz,u), i.e. V(—z,u) = V(z,u)C. The proof of Eq. (I63l) follows the line of reasoning used to 
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prove the factorization formula (1351) for the modular double in Sect. [H It relies on the properties ([5^ . 
(j55p of the elliptic gamma function. 

In order to elucidate the formula (j63p we indicate the matrix factors that are involved in Eq. (|63p 
at n = 1 and n = 2. Diagonal matrices 

Dy) = -^diagCe'^^ , 


D(2) = 


9i{2z - 2)ei{2z)9i{2z + 2) \ 


dmg(9i{2z-2)e‘^^^,-^^^9i{2z),9ii2z + 2)e-^^A 


9i{2rj) 


Matrices 


+ i) 


-03 {z - u) -03 {z + u) \ 
04 {z — u) 04 (z + tt) y 


/ _03 {z -u)9j{z - u + 2ri) 9s {z - u) 03 {z + u) 9s {z + u)9s{z + u - 2r/) \ 

^( 2 )+ i) = 9ys {z-u)9j^{z-u + 2r]) 0{3 {z - u) 04 } {z + u) 0^3 {z + u) 04 } (z + M - 2r/) 

\ 04 (z — n) 04 (z — n + 27 ]) 04 (z — u) 04 (z + u) 04 (z + n) 04 (z + u — 2 ?]) / 


The curly brackets in the second line of the previous formula denote symmetrization with respect to the 
theta function indices. Let us recall that at n = 1 the bases e and f, Eq. (1601) . are identical and the 
restricted R-operator coincides with the quantum elliptic Lax operator. The factorization of the elliptic 
L-operator appeared before in [25] . 
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